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We study numerically two versions of the monopole creation operators proposed by 
Frohlich and Marchetti. The disadvantage of the old version of the monopole creation 
operator is due to visibility of the Dirac string entering the definition of the creation 
operator in the theories with coexisting electric and magnetic charges. This problem does 
not exist for the new creation operator which is rather complicated. Using the Abelian 
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Higgs model with a compact gauge field we show that both definitions of the monopole 
creation operator can serve as order parameters for the connnement-deconnnement phase 
transition. The value of the monopole condensate for the old version depends on the 
length of Dirac string. However, as soon as the length is fixed the old operator certainly 
discriminates between the phases with condensed and non-condensed monopoles. 



: 1 Introduction 



The order parameters are useful tools for investigation of phase transitions. Up to now there 
is no good definition of the order parameter for the temperature phase transition in QCD with 
dynamical quarks. The traditional order parameters like the expectation value of the Polyakov 
line and the tension of the chromoelectric string work well in the quenched case (no dynamical 
quarks) but they fail to distinguish between the phases of the full QCD. Indeed, even at zero 
temperature the string spanned on a quark and an anti-quark can be broken by sea quarks 
making both the Polyakov line and the string tension useless candidates for the order parameter. 
Moreover, the value of the quark condensate (ty q ty q ) is strict order parameter only for massless 
dynamical quarks. 

Below we discuss the quantities which may serve as the order parameters for full QCD if 
the monopole (or, "dual superconductor") confinement mechanism |lj is valid. In this picture 
the monopoles defined with the help of the Abelian projection B are supposed to be condensed 



in the confinement phase. The monopole condensate causes a dual analogue of the Abrikosov 
vortex to be formed between quarks and anti-quarks. As a result the quarks and anti-quarks 
are confined into colorless states. In the deconfinement phase the monopoles are not condensed 
and quarks are not confined. Thus the natural confinement-deconfinement order parameter is 
the value of the monopole condensate, which should be nonzero in the confinement phase and 
zero in the deconfinement phase. 

There are two (formal) difficulties in defining of the monopole condensate. At first, the 
monopole condensate is defined as an expectation value of the monopole field. However, the 
monopoles are the topological defects in the compact Abelian gauge field and, as a result, 
the immediate output of the lattice simulations can not provide us with values the monopole 
fields. This difficulty can be overcame by noticing that the lattice output comes in a form of 
information about the monopole trajectories. Then one can apply a known procedure which 
allows us to rewrite the path integral over monopole trajectories as an integral over the monopole 
fields, and get the monopole condensate. 

The second difficulty is that the expectation value of the (e.g., scalar) field </> should al- 
ways be zero regardless whether this field is condensed or not. The reason is very simple: 
the path integral includes the integration over all possible gauges while the charged field is a 
gauge dependent quantity. These two problems were solved in Ref. ||, where the gauge in- 
variant monopole creation operator for compact QED (cQED) was explicitly constructed. The 
numerical calculations in the lattice cQED J|] and in the Maximal Abelian projection of the 
lattice SU(2) gluodynamics || confirm that the expectation value of this operator is an order 
parameter for the confinement-deconfinement phase transition. Other types of the monopole 
creation operators were investigated in Refs. |J. 

However, it noticed recently in Ref. || that the "old" monopole creation operator || de- 
pends on the shape of the Dirac string in the models with electrically charged dynamical fields. 
Exactly the same situation appears in the Abelian projection of QCD: the off-diagonal glu- 
ons become electrically charged dynamical fields while the diagonal gluons turns into compact 
Abelian gauge fields containing monopole singularities. Below we figure out whether the de- 
pendence on the Dirac string in the "old" version of the monopole creation operator is crucial 
for the usage of the expectation value of this operator as an order parameter. Another ( "new" ) 
monopole creation operator was suggested in Ref. [§]. This operator does not depend on the 
shape of the Dirac string even in the presence of the dynamical electric charges. However, the 
numerical calculations [0 show that the construction of this operator is rather complicated and 
the simulations are time consuming. 

The structure of the this paper is as follows. We explicitly describe the "new" and the "old" 
monopole creation operators in Sections 2 and 3, respectively. The results of the numerical 
calculations of these creation operators in the compact Abelian Higgs model are presented 
in Section 4. We demonstrate that both these operators can be used to detect confinement- 
deconfinement phase transitions in the theories with charged matter. Our conclusions are 
presented in Section 5. 
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2 "Old" monopole creation operator. 

The gauge invariant creation operator $ was suggested by Dirac |J: 

$ = <f>(x) exp |z J E k {x - y)A k {y) d 3 y^ , (1) 

here <p(x) and A k {x) are the electrically charged field and the gauge potential, which transform 
under the gauge transformations as 

<f>(x) -> 0(x)e iQ ^ , A k (x) -> A fc (ar) + <9 fc a(x) . (2) 

The Coulomb field, E k (x), satisfies the equation: 

d k E k = 5®(x). (3) 

It is easy to see that the operator $, eq. (P, is invariant under the gauge transformations (fj). 

The Frohlich-Marchetti construction H of the monopole creation operator in cQED is based 
on eq. ([!]). At first step the partition function of cQED is transformed to a dual representation. 
For the general form of the cQED action it can be shown || that the dual theory is an Abelian 
Higgs model (AHM) in the limit when both the Higgs boson mass and the gauge boson mass are 
infinite. In this theory the Higgs field, <f) x , corresponds to the monopole defect in the original 
cQED. The gauge field *B is dual to the original gauge field 9. Thus the gauge invariant 
creation operator (P for the AHM model corresponds to the monopole creation operator in the 
original cQED. The explicit expression for this operator on the lattice is (cf. eq.(Q)): 

$ mon,old = ^ e i(*B*H x ) ^ ^ 

where *H X is the Coulomb field of the monopole, 5*H X = *5 X , and *5 X is the discrete 5-function 
defined on the dual lattice. Here and below we will use the differential form notations on 
the lattice: (a, b) = Y2 C a cb c is the scalar product of the forms a and b defined on the c-sells; 
(a, a) = | \a\ \ 2 is the norm of the form a; d is the forward derivative (an analog of the gradient); 
5 is the backward derivative (an analog of the divergence) and *-operation transfers a form to 
the dual lattice. For a description of the language of the differential forms on the lattice see, 
e.g. 
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Performing the inverse duality transformation for the expectation value of the creation 
operator (|]) we get the expectation value, ($ mon ), of this operator in cQED: 

($ mon ' old ) = - y ve ex P {-s(de + w)} , z=j ve exp{-s(de)}, (5) 

where d9 is the plaquette angle, and S is the periodic lattice action, S^d^ + 27m) = S(d9), 
n G Z. The form W = 2n5A~~ 1 *(H x — uj x ) depends on the Dirac string *u x defined on the 
dual lattice. The Dirac strings start and end on the monopoles and anti-monopoles, 5*u x = 
*5 X . The operator A = d5 + 5d is the lattice Laplacian. The numerical investigation of this 
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creation operator in cQED shows [||] that it can be used as the confinement-deconfinement 
order parameter. 

The operator (Q) is well defined for the theories without dynamical matter fields. However, if 
an electrically charged matter is added, then the creation operator (P depends on the position 
of the Dirac string. To see this fact let us consider the compact Abelian Higgs model with the 
Villain form of the action: 

Zahm = I T>9 I V<p expj-/3||d# + 27m|| 2 - 7 ||dy? + g# + 27rZ|| 2 }. (6) 

Here 9 is the compact Abelian gauge field and <p is the phase of the dynamical Higgs field. The 
integer q is the charge of the Higgs field. For the sake of simplicity we consider the London 
limit (the Higgs mass is infinitely large while the Higgs condensate is fixed). 

Let us perform the Berezinsky-Kosterlitz-Thouless (BKT) transformation [IT] with respect 
to the compact gauge field 9: 

d9 + 2vm = dA + 2Tr6A~ 1 j , with A = 9 + 2n5A- 1 m[j] + 2vr£; . (7) 

Here A is the non-compact gauge field, *m[j] is a surface on the dual lattice spanned on the 
monopole current *j (5*m[j] = *j) and k is the integer-valued vector formQ. We substitute 
eqs.([?D in eq.@ and shift of the integer variable, I — > I + qk. 

Next we perform the BKT transformation with respect to the compact scalar field ip: 

dp + 2irl = d$ + 27r6A~V , with & = if + 2tt(5A" 1 s[(t] + 2np . (8) 

Here $ is the non-compact scalar field, *s[a] is a 3D hyper-surface on the dual lattice spanned 
on the closed surface *a (S*s[a] = *a) and p is the integer-valued scalar form. 

Substituting eqs. ([7],j8p into the partition function (|6]) and integrating the fields A and ip we 
get the representation of the compact AHM in terms of the monopole currents *j and world 
sheets of Abrikosov strings V ("the BKT-representation"): 

Zahm oc Z bkt = £ £ exp{ - 4tt 2 /3 (j, A + m2 j) ~ 4tt 2 7 (a,, A + m 2 ^) } ' ( 9 ) 

*j€Z(*c 3 ) * CT;/ eZ(*c 2 ) 
S*j=0 S*<jj=q*j 

the dual surface variable Vj = *cr + q*m[j] is spanned g-times on the monopole current j, 
5*<jj = q*j, since the flux of the magnetic monopole having an unit magnetic charge can be 
taken out by q strings carrying the elementary flux 27r/g. The mass of the gauge boson 9 is 

The BKT-representation (^) of the AHM partition function ([]) can also be transformed 
into the dual representation using simple Gaussian integrations. We use two dual compact 
fields *B (vector field) and *£ (scalar field) in order to represent the closeness properties of the 

1 A detailed description of the duality and BKT transformations in terms of the differential forms on the 
lattice can be found, e.g., in Ref. [jT~o|| . 
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currents *<x, and *j, respectively. We also introduce two dual non-compact fields, *F (vector 
field) and *G (rank-2 tensor field), in order to get a linear dependence on, correspondingly, the 
currents *Oj and *j under the exponential function: 

/oo /*oo pit Pit 

V*F / V*G / V*B / V*i Yl 
J- 7 * *jez(* C3 )*a J ez(* C2 ) 

exp {*G, (A + m 2 )*G) - * 7 (*F, (A + m 2 )*F) 

+i{*F, V,-) + i(*G, *j) + i(*B, 5*a 3 - q*j) - i(% 5*j) } , (10) 

where 

* (3= m^' * 7= i6^' 

Note that in this representation the integer variables *<jj and *j are no more restricted by the 
closeness relations. Therefore we can use the Poisson summation formula with respect to these 
variables and integrate out the fields *F and *G. Finally, we obtain the dual field representation 
of the partition function (|[): 

V*B / V*£ ^ 

exp{-*/3 (d*B + 2n*u, (A + m 2 ) (d*B + 2n*u)) (12) 
-*7 (d*^ + <f # + 2tt*i;, (A + m 2 ) (d*£ + q*B + 2vr*t;)) } , 

where *u and *t> are the integer valued forms defined on the plaquettes and links of the dual 
lattice, respectively. Clearly, this is the dual Abelian Higgs model with the modified action. 
The gauge field *B is compact and the radial variable of the Higgs field is frozen. The model 
is in the London limit and the dynamical scalar variable is the phase of the Higgs field *£. 

Thus in the presence of the dynamical matter the dual gauge field *B becomes compact^. 
The compactness of the dual gauge field implies that it transforms under the gauge transfor- 
mations as: 

*B -> *B + d*a + 2n*k , (13) 

where the integer valued field k is chosen in such a way that *B e (— tt, it}. 

One can easily check that the operator (|j) is not invariant under these gauge transforma- 
tions: 



2 Another way to establish this fact is to realize that the pure compact gauge model is dual to the non- 
compact U(l) with matter fields (referred above as the (dual) Abelian Higgs model). Reading this relation 
backwards one can conclude that the presence of the matter field leads to the compactification of the dual 
gauge field *B. 
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3 "New" monopole creation operator. 

The invariance of the operator @) under the gauge transformations (|13"D can be achieved if and 
only if the function *H X is an integer- valued form. Thus, if we take into account the Maxwell 
equation, 5*H X = *5 X , we find that *H X should be a string attached to the monopole ("the 
Mandelstam string"): *H X — > *j x , *j x G Z, 5*j x = *5 X . The string must belong to the three- 
dimensional time-slice. However, one can show || that for a fixed string position the operator 
$ creates a state with an infinite energy. This difficulty may be bypassed || by summation 
over all possible positions of the Mandelstam strings with a measure /i(*j): 

jxEZ 
S*jx = *S x 

The summation over the strings provides an entropy factor which cancels the energy suppres- 
sion. An example of a "reasonable" measure (i(j x ) is ||: 

n(*j x ) = ex p{~^ir^n 2 } • ( l6 ) 

If the Higgs field <fi is g-charged (q G Z), the summation in eq. fllSl) should be taken over q 
different strings with the unit flux. Measure (|16D corresponds to the dual formulation of the 
3D XY-model with the Villain action: 

S XY (x,r) = |||d X - 2nB + 2nr\\ 2 . (17) 

The Mandelstam strings correspond to the ordinary vortices in the XY-model (|Hj). 
The two point correlation function in the XY-model (0) is given by 



t^y Jxyl ' 

Jxy 

where 5*j xy = *5 X — *5 y . For sufficiently large coupling k and sufficiently small field B we get 

(e iXx e~ iXv ) -»• const , (19) 

as | a; — y\ — > oo. So that the correlation function (|18|) might yield an appropriate measure of 
the Mandelstam strings. 

Performing the inverse duality transformation of the creation operator (|15D in the dual 
representation of the compact AHM (12) we get the expectation value of this operator in terms 
of the compact gauge field 9 in the original representation (H): 



($ 



moil, new 



) = | E /^^exp{-i-||d*j|| 2 -5(d^ + ^j)}- (20) 



*jez(*c 3 ) 
5*j=*8x 



The current j = *( 3 )*( 4 )j is defined via a double duality operation which is first applied in the 
3D time slice and then in the full 4D space. 

We thus defined "old" (EI) and "new" (|20D monopole creation operators. 
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4 Numerical results 



Below we present results of the numerical simulation of the old and the new monopole creation 
operators. We study them in the Abelian Higgs model with compact gauge field and infinitely 
deep potential for the Higgs field corresponding to the London limit. This is the simplest model 
containing both the monopoles and the electrically charged fields. We used the Wilson form of 
the action which is more suitable for the numerical simulations than the Villain action used for 
the analytical calculations above: 

Zj HM = I V9 I v^e Pcosm+lcos{A ^ qe) . (21) 

J —IT J —TV 

The expectation values of the old and new monopole creation operators in the Wilson form are 
given by 

( $ mon,oid^ = 1 f D o exp {[] cos (d6 + 27i5A- 1 *(H x -u x ))+ 1 cos(q6)}, (22) 

Z J-TT 

($ mon ' ncw ) = i [ ^exp{--^||d*j|| 2 + /?cos (d6 + — j) + 7 cos(g#)}, (23) 

8*j=*S x 

respectively. We have fixed the unitary gauge, therefore the phase of the Higgs field is eaten 
up by the corresponding gauge transformation. 

The value of the monopole order parameter, ((f)) , corresponds to the minimum of the (effec- 
tive constraint) potential on the monopole field. This potential can be estimated as follows: 

V eS ($) = - In (s($- $ mon > new /oidj ^ . (24) 

We are studying the model with q = 7. The compact Abelian Higgs model with multiple 
charged Higgs fields is known to have a nontrivial phase structure containing the Coulomb, 
Higgs and Confinement phases |T2[. In this paper we concentrate on the phase transition 



between the confinement and the Coulomb phases. 

We simulated the 4D Abelian Higgs model with anti-periodic boundary conditions in the 
3D space (the single monopole charge can not exist in the finite volume with periodic boundary 
conditions). For the new creation operator we perform our calculations on the 4 4 , 6 4 , 8 4 lattices 
and the coupling constant 7 = 0.3. For each configuration of 4D fields we simulated 3D 
model to get the ensemble of the Mandelstam strings with the weight fi{j x )- We generated 
60 statistically independent 4D field configurations, and for each of these configurations we 
generated 40 configurations of 3D Mandelstam strings. 

According to eq. (]17|) the weight function (|16|) corresponds to the 3D XY-model with the 
Villain action. The XY model has the phase transition at k c (B = 0) ~ 0.32 ||13|| . Our numerical 
observations show that in presence of the external field B the critical coupling constant gets 
shifted, «f ^0.42. 

Typical configurations of the Mandelstam strings corresponding to the phase with the con- 
densate of these strings (disordered phase, large k) and to the phase without the string con- 
densate (ordered phase, small k) are shown in Figures |l|(a,b). 
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(b) 



Figure 1: The typical Mandelstam strings in the auxiliary 3D model in (a) the ordered phase 
(no string condensate, k = 0.3) and (b) disordered phase (nonzero string condensate, k = 0.5). 



One can expect that the operator fll5|) plays the role of the order parameter in the phase, 
where the Mandelstam strings are condensed (k > k c ). In Figures we present the effective 
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Figure 2: The effective monopole potential 
phases. 



4j) in (a) confinement and (b) deconfinement 



potential (|23J) in the confinement {(3 = 0.85) and deconfinement (f3 = 1.05) phases. The 
potential is shown for two values of the 3D coupling constants k > k c corresponding to high 
densities of the Mandelstam strings. In the confinement phase, Figure 0(a), the potential V($) 
has a Higgs form signaling the monopole condensation. According to our numerical observations 
this statement does not depend on the lattice volume. In the deconfinement phase, Figure §(b), 
the potential has minimum at $ = which indicates the absence of the monopole condensate. 

For small values of the 3D coupling constant k (in the phase where Mandelstam strings 
j x are not condensed), we found (Figure |3|) that the potential V(<E>) has the same behavior 
for the both phases of 4D model. Thus the operator (O) serves as the order parameter for 
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Figure 3: The effective monopole potential ( P4|) in the low— k region of the 3D model. 

the deconfinement phase transition, if Mandelstam strings are condensed, i.e. the coupling 
constant k should be larger than the critical value k c (B). 

Now let us show that the expectation value of the old monopole creation operator (||) 
behaves as an order parameter for the Dirac string with a fixed length. We generate the 
compact Abelian Higgs model for 7 = 0.3, q — 1 and couplings (3 = 0.6 (confinement phase) 
and (3 = 1.2 (deconfinement phase). 

We use Dirac strings of the given length in proportion 1 : 1.5 on the 4 4 , 6 4 , 8 4 , 10 4 , 12 4 , 16 4 
lattices. The form of the Dirac strings is shown in Figure £|. For each value of (3 and for 
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Figure 4: Two configurations of the Dirac strings used in our simulations on the lattice with 
anti-periodic boundary conditions. 
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the fixed form of the Dirac string we took average over 3000 values of the monopole creation 
operator. The dependence of the minimum of the effective potential, minV($), on lattice of 
the size L for the fixed length of Dirac string i is shown in Figure || We fitted the data for 



min V(<£) 




o — ^ — - — ■ ■ ■ ■ ^ 

0.05 0.1 0.15 0.2 0.25 



1/L 

Figure 5: The minimum of effective monopole potential ([24] ) as a function of the inverse lattice 
size 1/L for the Dirac strings with various lengths i. 

^min by the formula $ m i n = aL b + where a, b are fitting parameters. It occurs that b = — 1 
within statistical errors. The resulting values of are collected in Table 1. We obtain that 
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-0.03(6) 
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-0.07(5) 



Table 1: The minimum of the monopole potential, $^in5 based on the old operator vs. the 
gauge coupling (3 and the Mandelstam string length £. 

^min vanishes at the point of the phase transition. Thus the operator (|15"D serves as the order 
parameter for the deconfinement phase transition provided the length of Dirac string is fixed. 
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5 Conclusions 



The new monopole creation operator proposed in Ref. || can be used as a test of the monopole 
condensation in the theories with electrically charged matter fields. Our calculations indicate 
that the operator should be defined in the phase where the Mandelstam strings are condensed. 
The minimum of the effective potential corresponds to the value of the monopole field which is 
zero in deconfinement phase and nonzero in the confinement phase. 

In the infinite volume limit the potential corresponding to the old version of the monopole 
creation operator shows the same features as the potential calculated with the use of the new 
operator. The shape of the old effective potential depends on the length of the Dirac string. This 
fact indicates that the Dirac strings with different shapes provide different monopole creation 
operators and all of them can serve as order parameters for the confinement-deconfinement 
phase transition. 
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